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Abstract — In this paper we will explore fundamental con- 
straints on the evolution of certain symplectic subvolumes 
possessed by any Hamiltonian phase space. This research has 
direct application to optimal control and control of conservative 
mechanical systems. We relate geometric invariants of sym- 
plectic topology to computations that can easily be carried out 
with the state transition matrix of the flow map. We will show 
how certain symplectic subvolumes have a minimal obtainable 
volume; further if the subvolume dimension equals the phase 
space dimension, this constraint reduces to Liouville's Theorem. 
Finally we present a preferred basis that, for a given canonical 
transformation, has certain minimality properties with regards 
to the local volume expansion of phase space. 



I. Introduction 



A. Overview 



The traditional approach for studying the dynamics and 
control of mechanical systems is to focus on individual 
trajectories and states in order to determine where they will 
go and where they can be forced to go. In reality, however, 
system states are never precisely known and can only be 
determined to exist within some set of finite volume in the 
dynamical system's phase space. By treating such systems 
as a sum of individual trajectories, one loses the geomet- 
rical insight and deeper results offered by more wholistic 
approaches. 

In this paper, we will be concerned with understanding 
fundamental constraints on the evolution of compact 2k- 
dimensional symplectic sets that evolve along a nominal 
trajectory of the system. Different symplectic constraints 
arise on such sets, including conservation of the signed 
2 fc- volume projections on the coupled symplectic planes as 
well as the constraints implied from Gromov's Nonsqueezing 
Theorem (see Scheeres et al [13] for a discussion of these 
constraints in relation to orbit uncertainty evolution). We 
will further present an additional constraint for a minimal 
obtainable volume that exists on certain classes of 2k- 
dimensional symplectic sets and show how such a constraint 
leads to the local collapse of phase space along solution 
curves in Hamiltonian phase space. This collapse of phase 
space is fundamentally linked to the expansion of symplectic 
subvolumes. Finally, for any fixed final time, no matter how 
large, we will produce a distinguished orthogonal symplectic 
basis that resists collapse. The basis may collapse as time 
evolves, but will return to being orthogonal at the final time. 
The uncertainty of any 2fc-dimensional distributions initially 
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parallel to the symplectic planes of this basis, even though it 
may increase dramatically during the course of its evolution, 
will always return to its initial uncertainty at the final time. 

Since the resulting equations produced by applying Pon- 
tryagin's Maximum Principle to optimal control problems 
are Hamiltonian, the results we discuss here should provide 
geometric insight to the evolution and control of uncertainty 
distributions in such systems. This theory provides funda- 
mental limits on dynamical orbits, and hence if one provides 
a control it provides limits on the accuracy of the control 
in the face of uncertainty. It also provides limits on uncer- 
tainty propagation in optimal control systems. Moreover, the 
preferred minimal uncertainty basis we produce should have 
numerous benefits to the design of fixed finite time optimal 
control laws where precise state information is unknown. 

B. Outline 

In Section 2, we introduce Hamiltonian systems and 
the state transition matrix (STM). We show how classical 
identities on the Lagrange and Poisson brackets relate to 
constraints on the STM. Specifically, for any symplectic 
column of the STM $, the sum of the 2x2 symplectic 
subdeterminants must add up to unity. 

In Section 3, we look in depth at surfaces that can be 
explicitly parameterized by one of their symplectic planes. 
We derive area expansion factors from the parameterization 
plane to the surface, its image under the Hamiltonian phase 
flow, and the symplectic projections of its image. If the state 
of the system is somewhere on the initial surface (with equal 
a priori probability), we interpret these various expansion 
factors as a probability map that leads one to understand 
where the particle, after applying the Hamiltonian phase flow, 
is likely to be found. 

In Section 4, we present a nonsqueezing-like property for 
2fc-volumes which is closely related to Wirtinger's Inequality. 
We show how this leads to the fact that subvolume expan- 
sions in the differential neighborhood of the Hamiltonian 
flow leads to a collapsing property in systems which exhibit 
chaos. 

In Section 5 we will discuss how these constraints on 
subvolume expansions, when considered with Liouville's 
Theorem, leads to the local collapse of the phase space 
around nominal trajectories. Interestingly, we will also show, 
that given a canonical transformation, there exists a preferred 
basis that resists collapse. In particular, the volume of a 
2k subvolume chosen to be initially parallel to k of the 
symplectic planes will return to its initial value at this fixed 
final time. 



II. Hamiltonian Systems 

A. The Classical Approach 

Hamilton 's Equations: In an N degree of freedom Hamil- 
tonian System, one has a 2Af -dimensional phase space 
spanned by generalized coordinates {qi}i and their 
conjugate momenta {pi}^ ■ The conjugate pairings of 
coordinates and momenta, {qi,Pi) form what are known 
as symplectic pairs. The dynamical equations of motion 
are deriveable from a Hamiltonian function H{q,p) and 
Hamilton's equations: 

qi = ^ and Pi ^ (1) 

opt dq., 

Solutions curves of the system ([U are called the Hamiltonian 
phase flow, and are denoted (l)t{q,p). 

Matrix Formalism: Let x = (pi, qi,P2, 92, • ■ -Pn, qn) be 
the phase space position of the system, ordered by symplectic 
pairs. Define the matrix 

■ -1 
1 



J2 



The symplectic matrix J is defined as the 2N x 2N block- 
diagonal matrix with J2's down the main diagonal. Grouping 
the coordinates in symplectic pairs will be useful to us later 
as we will be looking at various "symplectic columns" of 
the STM. Given the above definitions, Hamilton's equations 
can be cast into the following matrix form: 



at ox 



(2) 



where x = <j>t{xo). The right hand side is the so-called 
symplectic gradient of the Hamiltonian function. 
Lagrange Brackets: Given a transformation 

Qt^Qt{q,p) and Pi = Pi{q,p) 

we may introduce the Lagrange bracket expression for the 
two variables (which can take on any of the values 

qi, . . . ,9„,pi, . . . ,p„): 

[u,v] = 2^\—— — (3) 



V du dv 



du dv 



The exactness conditions required for a canonical trans- 
formation can then be cast into the following equivalent 
conditions: 



, 9fc] = [pj,Pk] = [Pj,(lk] 



(4) 



The Hamiltonian phase flow (f>t{q,p) is a continuous one 
parameter family of canonical transformations. 

Poisson Brackets: Alternatively, one can define the Pois- 
son bracket as: 



N 



du dv du dv 
dpi dqi dqi dpi 



(5) 



where (u, v) can now be any of the variables 
Qi, . . . ,Qn, Pi, ■ ■ ■ , Pn- The sufficient conditions for 
a canonical transformation can also be written as follows, 
in terms of the Poisson bracket: 



B. The Geometric Approach 

Symplectic Manifolds: A symplectic structure on an even- 
dimensional manifold M is a closed nondegenerate differen- 
tial two-form uj on M: 



and 



The form u; is called the symplectic form and the pair (M, oj) 
is called a symplectic manifold. 

On any symplectic manifold, there exists a vector space 
isomorphism between its cotangent and tangent bundles. At 
X e M, we have 

4 : T;M T,M 

defined by the symplectic form and the following relation. 
A vector ^ e T.j.,M is mapped to the one-form which 
acts on a vector rj E T^M as follows: I^^iOiv) — ^{ViO- 

Hamiltonian Flows: Let H he a function H : M ^ M. 
which we will call the Hamiltonian. The associated Hamil- 
tonian vector field on M is defined by IdH. The flow 
generated by the vector field IdH is the Hamilton phase 
flow (pf If M = with the standard symplectic form 
ujQ = dpi A dqi, we recover Hamilton's equations ([T]i. 

A transformation : A/ M is considered symplectic or 
canonical if it preserves the symplectic form, i.e., (ffuj = uj. 
The Hamiltonian phase flow (f>t is a one parameter family of 
canonical transformations. 

Integral Invariants: A differential /c-form a is an integral 
invariant of the map (j> if the integrals of a on any fc-chain 
a is preserved as follows: 



The symplectic form uj is an integral invariant of the Hamil- 
tonian flow. 

To gain a physical intuition for what the integral invariant 
J UJ represents, consider now a closed parametrized surface 

-!/;(cr 
(T C 



m 

n>2 _ 



^" ^ (p, with a parametrization given by ■0 : 

■0 : [u^v) — > {q{u,v),p{u,v)). Then 



dp A dq = 



i=l 



d{pi,qi) 
d(u, v) 



du dv 



Ilcr %luv) dudv represents the oriented area of the pro- 
jection of the surface 0(ct) on the i-th symplectic plane. 
By considering 0((t) as an initial surface and applying the 
Hamiltonian phase flow, mapping the surface to 04(?/'(cr)); 
we recognize, as a physical interpretation of the preservation 
of the symplectic form under canonical mappings, that the 
sum of the oriented areas of the projections onto the TV 
symplectic planes is preserved. 

C. The State Transition Matrix 

Definition: If : Af ^ M,(j}{p,q) = {P,Q) is a 
canonical transformation, its differential 



{Qj, Qk} = {Pj,Pk} = {Pj,Qk} = 6jk 



is, when represented in matrix form, known as the state 
transition matrix (STM) (f>, a terminology adopted from 
Linear Systems Theory. Supposing (j>t sends (f)t : xq ^ 
4>{xo), the STM <f> maps initial deviations in the initial 
conditions $ : Sxq ^ ■ 6xq to its final state, so that 
(f)t {xq + Sxq) w X + $ • 6x0, to first order. Here x = {p, q) G 
M. 

Dynamics: If (j>t is the Hamiltonian phase flow, we have: 
If we perturb the initial conditions to xq 



5xo, we find: 



(6) 



-^4>t{xQ + 5xq) = J • ^^{x + 5x) 
By expanding this in a Taylor Series, one sees: 

— $ — J ■ • <!> 

dt 

This defines a system of AN"^ differential equations that 
can be integrated numerically, simultaneously along with the 
nominal solution curve 0t(a;o). 

Relation to Lagrange and Poisson Brackets: We will relate 
the Lagrange and Poisson Brackets to determinants of vari- 
ous submatrices of the STM. We will arrange the coordinates 
in a symplectic order, so that x ~ {pi,qi, ■ ■ ■ ,Pn, 9Ar) and 
X = (Pi, Qi, . . . , P„, Qat). In this fashion, the STM $ is 
thought of as 

ox 

Define the following subdeterminants: 



Mij — det 



dPj dPj ~ 



Hence, Mij is the subdeterminant of the intersection of 
the ith symplectic row of the STM with its jth symplectic 
column. 

It is easy to see the Lagrange and Poisson brackets are 
related to these subdeterminants as follows: 



N 

i=l 
N 



(7) 



(8) 



III. Symplectic Surfaces 

A. Surface Classifications 

We will begin by making the following fairly natural 
definitions. 

Definition 1: A globally symplectic surface is a two- 
dimensional submanifold of the phase space R^'i which 
admits a 1-to-l mapping to at least one of the symplectic 
planes via the projection operator, i.e. it is a surface which 
can be parameterized in explicit form by one of its symplectic 
coordinate pairs. 

This characteristic is not an invariant one. It is possible, 
for example, for a lamina parallel to a symplectic plane to 



fold under some symplectic map, so that its image under the 
map is not 1-to-l with any symplectic plane. 

Another surface type we will consider is the following: 
Definition 2: A parametrically symplectic surface (or 
parasymplectic surface, for short) is a 2-dimensional sub- 
manifold of R-^n that admits a parameterization that is a 
symplectic one, i.e. one with a parameterization map that 
is canonical. 

The parasymplecticity of a surface is an invariant char- 
acteristic. Let (7 be a lamina on the symplectic plane (u, v) 
which is the parameterization of the surface (f){<j) C M^", 
where is the parameterization map. Let be a symplecto- 
morphism, which exists if (t){a) is a parasymplectic surface. 
Let ij) : R2" R2" be an arbitrary symplectomorphism 
which takes i/j : (t){a) — > ■ip{(j){a)). Then ■i/;((/)(cr)) is 
parasymplectic, with symplectic parameterization ?/;o0 : cr 

We will consider tt^ : R^" _^ jja symplectic pro- 

jection operator, so that ni{{pi,qi, ...,pn, qw}) = {Pi, qi)- 

B. Area Expansion Factors 

Notation: We define the 2n x 2 matrix 11^ as: 

= [ O2 O2 • • • O2 I2 O2 • • • 02]^ 

where the I2 appears in the K-th symplectic row. For any 
2n X 2n matrix A, the product A ■ 11^ is the K-th symplectic 
column of A; 11^ • A is the K-th symplectic row of A; and 
• A ■ Ii\ is the 2x2 intersection of the K-th symplectic 
row with the A-th symplectic column. 

In this notation, the STM subdeterminant Mij, defined 
previously, can be expressed as: 



det(nf • $ • Hj) 



Globally Symplectic Surfaces: We will consider now a 
surface r which is globally symplectic with respect to the j- 
th symplectic plane; i.e., the projection map TTj : t C R^" ^ 
7rj(T) C R2 is one-to-one. We can parameterize t by its 
symplectic shadow on the j-th symplectic plane. Now let the 
surface r be mapped into the future by the Hamiltonian flow 
(f)t ■■ . . . ,_pw,gw) {Pi,Qi,. . . ,Pn,Qn)- We will 

now consider the projection of (f)t (r) onto the i-th symplectic 
plane. For our analysis we will consider a differential area 
element dr of t (see Fig. [T]). 

The surface r is described by the 

Pj and V — qj 
■ ■ ,PN{u,v),qN{u,v)). 
as follows: 



rametenzation u 
{pi{u,v),qi{u,v),. 
define the matrix L 



,u,v, 

dxj 



pa- 
by 
We 



L 



dpi 
du 

dqi 
du 



du 
aqN 
du 



dv 



dpN 
dv 

aqN 
dv 




Fig. 1. Area Expansion Factors 



where the j-th symplectic row is equated to I2. L is the 
Jacobian matrix of the parameterization map that takes 
TTj{dT) — > dr. The metric for the surface r in terms of 
variations in the u — v plane is given by: 



dx dx. dx dx 



dv du dv dv 



where x = {pi{u,v), . . . ,qN{u,v)), as usual. Hence the 
metric determinant is: 

g = det{g,j) = ©(£) = det(L^ • L) 

where is the Gram determinant of the matrix L, defined 
by this equation. 

The surface area of r is thus given by: 



S.A. 



y/gdudv 



A simple application of this result gives us the total 
physical area of the area element dr in terms of the area 
of its projection: 

A{dT) 



A{n,{dT)) 

Similarly we have: 

A{7T,{dT)) - 



where $ is the STM associated with (pt and xq- 

The area of the projection TTi{(f>t{dT)) is given by the 
Jacobian: 

djPi, Qi) _ / dl^dQ^ _ dQ^dPA dx^ dxx 

d{u,v) ^^^^\dxK,dx\ Oxk dxx J du dv 

which can be represented more concisely as: 

Ai7:,{MdT))) 



A{nj{dT)) 



= dct(nf • $ • L) 



Pamsymplectic Surfaces: We will now consider the case 
where t is parallel to the j-th symplectic plane. In this case, 
its parameterization map is a symplectic one, and thus it is 
a parasymplectic surface. All of the above results hold, but 
the matrix L reduces to the simpler form L = Hj, which 
gives us the following: 

A{dT) =A{TTj{dT)) 



AiMdr)) 

A{TTj{dT)) 

and, most notably 

AiMMdr))) 



©(<!>• n,) 



Ain^idr)) 



det(nf • $ • Hj) = Mi. 



Preservation of the sum of the oriented symplectic area 
projections thus gives us the following constraint on the 
STM: 



N 



= [pj,qj] = 1 



(9) 



Application to Orbit Uncertainty Distributions: Suppose 
we know that a system can be found anywhere on the surface 
T with equal a priori probability. The surface is now mapped 
into the future by the Hamiltonian flow (pt and we wish 
to determine where the particle is most likely to be on the 
Pi - Qi plane. 

We begin by discretizing the (u — Pj) — {v = qj) plane, 
each area element with area AuAt;. Summation will be 
assumed to be over each district. The probability that the 
particle is in dr is given by: 



V{dT) 



A{dT) y0(L)AwA?; _ 



A{t) J2 v/©(I) AuAt; E 76(1) 
We reiterate the area of TTi{ipt{dT)) is given by: 

Ain,{(j)t{dT))) « det(nf • $ • L)AuAv 
so that the area probability density at {Pi, Qi) is: 



I det(nf • $ • L)\AuAvJ2vW^ 

This approach may be helpful in asteroid tracking, where 
angular and angular rate information is precisely known, 
but there is initial uncertainty in the r, r distribution. This 
problem is treated in Milani, et al. [8]. 

IV. Integral Invariants 

In this section we will discuss the difference between two 
fundamental integral invariants defined for an arbitrary 2k- 
dimensional subvolume of our 2n-dimensional phase space, 
for k = l,...,n. The first is the well-known integral 
invariant of Poincare-Cartan. The second integral invariant 
is closely related, and is tantamount to a global version of 
Wirtinger's Inequality for lower dimensional subvolumes of 
phase space. 



A. Signed and Unsigned Integrals of Differential Forms 

We will begin our discussion on integral invariants with a 
brief discussion of the theory of integration of differential 
forms. Let S c M^" be a 2fc-dimensional submanifold, 
parameterized by </> : (a C M."^^) — > S. Let a be a 2fc-form 
on R^". Then we define 



a 



For some function f{xi, . . . , X2k)^ the puUback of a can be 
expressed in the following form: 



,2fe 



*a = f{xi, . . . , X2k)dx^ A - ■■ Adx 



where we take {xi,...,X2k) to be a basis of R^*^. Then 
the integral of <p*a over a reduces to the ordinary euclidean 
integral: 

a= 4'*a= I f{xi,...,X2k)dx'^ ■■■dx^'^ 

We would like to introduce a further definition as follows. 
We define the unsigned integral of a over S to be: 

/ \a\= \f{xu...,X2k)\dx'^ ■■■dx'^'' 

where / has been defined above. 

A more rigorous definition of the above integrals must 
involve a partition of unity, but simplicity has been choosen 
over rigor so as to illustrate the spirit of the definitions. 

B. The Integral Invariants of Poincare-Cartan 
Consider the standard symplectic form 



^PtAqi 



and its fc-th exterior product: 



1 k 

— CO 

k\ 



dpi^ A dqi^ A • • • A dpi^ A dqi^ 



l<'il <'"<U- <n 

Consider a set of 2k vectors (X^, . . . , X'^'') in R^" Then 

represents the sum of the oriented 2fc-volume projections of 
the parallelpiped spanned by X^, . . . , X'^'^ on the symplectic 
"2A;-planes." 

w*' is known as the integral invariant of Poincare-Cartan. 

Given an arbitrary 2fc-dimensional phase volume fl (in a 2n 
dimensional space) and the Hamiltonian phase flow (f>t, we 
have: 

so that the sum of the oriented 2A;-volume projections on 
each symplectic "2k plane" is conserved. 



C. The Wirtinger-Type Integral Invariants and Volume 
Identifying C" ~ R^", the symplectic form becomes 



= ^"^dziAdzi 



for any k € [l.N], consider the vectors Xi, . . . ,X2k S 
C" ~ M^". Wirtinger's Inequality states that the "2k" volume 
of the parallelpiped spanned by these 2k vectors is bounded 
by 



■ ,x2k)\<yohk{Xi 



, X2k) 



(10) 



We make the following two observations. First, it is clear 
that \u)^\ is an integral invariant of the Hamiltonian flow (pt, 
so that, given any 2A;-volume O, we have: 

k\ Jo ' ' k\ 



CO" 



Moreover, 



1 



< Vol2fe(r!) 



so that this integral invariant represents a minimum (2k) 
volume that the body fl may obtain. For the case k = n, 
the volume of O is a constant which equals this invariant 
quantity (Liouville's Theorem). 

D. Parasymplectic 2k-Volumes 

In direct analogy with our discussion of parasymplectic 
surfaces, we define parasymplectic 2fc-volumes as follows: 

Definition 3: A parasymplectic 2fc-volume, or paramet- 
rically symplectic 2fc-volume, is one that admits a parame- 
terization whose paramteterization map is a symplectic one. 

An example of a parasympletic volume is the following. 
Take any 2fc-dimensional volume that is parallel to k of 
the symplectic planes, i.e. a region defined by: 

^ = • • • ,Pk,qk,Ck+l,dk+l, ■ ■ ■ ,Cn,dn) 

where the variables pi,qi, . . . ,pk,qk vary over some region 
of M?'' and Ck+i, dk+i, ■ ■ ■ Cn, dn are constants. Now let 
be the Hamiltonian phase flow. The 2fc-phase volumes 
(l>t{^) C R^" are a one-parameter family of parasymplectic 
2fc-volumes. 

Theorem 1: [Volume Expansion of Parasymplectic 2k- 
Volumes] Let D. c R^'^ be the parameterization of a 

volume (t){Vl) C R^" in a symplectic phase space whose 
parameterization map is a symplectic one. Then 

V0l2fe(l^) < Y0\2k{(l>{^)) 

Corollary 1: The 2A;- volume of any parasymplectic 2k- 
volume is at least as large as the volume of its symplectic 
parameterization . 

Proof: To prove our theorem, we only need prove it 
for a differential volume element. The generalization foUows 
via a simple integration argument. 

Let (j) : R^'^ —^ R^" be a symplectic parameterization of a 
surface. Let 

• d d d d \ 

dui dvi ' ' duk ' dvk J 



be a basis of M^'^ and let 
d_ 

. dui 



and 



_d_ 

dvi 



be the push forwards of the basis vectors in the parameteri- 
zation space. Applying Wirtinger's Inequality ( fTOl i. we have 
that 



1 [A. A. 

kl \dui ' ' dvk 



— dui A • • • A dvk 



d d \ 



(9ui ' ' dvk ) 

- 1 < voi2fe(Xi,./.n) 

Hence any 2/c-dimensional volume measure must be nonde- 
creasing under such a map. ■ 

E. The Volume Expansion Factor 

In this section we provide a practical approach to de- 
termining the volume and the integral invariants of 2k 
sub volumes. We will consider the volume i7 C M^'"' to be the 
parameterization volume of a 2fc-volume in the symplectic 
space 



i.^", with parameterization map 



The Jacobian of the parameterization is the 2n x 2k matrix 
given by 

L = d(j) 

We will be interested in computing the total volume of fl, 
the sum of its oriented symplectic projections (i.e. integral 
invariant of Poincare-Cartan), and its minimum obtainable 
volume (i.e. the integral invariant of the Wirtinger type). 

In terms of the parameterization coordinates 
{ui,vi, . . . ,Uk,Vk) E M^'', the following metric is 
induced on the surface: 



g = . L 

We thus recognize the determinant of the metric g = det g 
as the Gram determinant of the Jacobian matrix L: 



g = e{L) = dct(i^ • L) 



so that 



Jn 



In practical terms, the Gramian of the Jacobian can be 
identified with the volume expansion factor: 

where V2k{dVl\ (/>) is the local 2fc-volume expansion factor of 
d^l under the mapping (j>. 



V. Local Collapse of Phase Space 

A. Volume Expansion and the Local Collapse of Phase Space 

The setting for this subsection will be the evolution 
of a differential neighborhood surrounding a Hamiltonian 
trajectory through phase space. Consider the Hamiltonian 
flow: 



(j)t : . . . ,Pn,qn) {Pl,Ql 



Now consider a differential 2n-"cube" D, situated at the initial 
point X = {pi,qi, . . .pn,qn), whose faces are parallel with 
the symplectic planes. Let T C 51 be a 2fc-dimensional subset 
that is parallel with k of the symplectic planes, and let T' C 
r2 be a 2n — 2fc dimensional subset that is parallel with the 
remaining n — k symplectic planes, such that 57 is a direct 
sum: 

O = T ® T' 

and, therefore 

Vol2„(r!) = V0l2fc(T) • V0l2„_2fc(T') 

The Hamiltonian flow now takes (j)t : x X along with its 
differential neighborhood. We define 

n = cbtm 
T = MV 

T = 0t(T') 
We now define the angle [3 via the relation: 

Yohnin) = Vol2fc(T) • Vol2„-2fc(T^)sin/3 

so that f3 is the angle between the subspaces T and T'. By 
Liouville's Theorem, we have: 



Vol2„(r!) - Vol2„(r!) 



so that: 



1 = i^2fe(T;0t)i^2n-2fe(T';(/)t)sin/3 

But by Theorem [T] we have: 

J^2fc(T;0t) > 1 
J^2«-2fe(T';0t) > 1 

We conclude that the greater the volume expansion of these 
lower dimensional differential "slices" T and T', the greater 
the inward collapse of their respective subspaces towards 
each other. In chaos theory, where Vol2fc((/)t(T)) is growing 
at an exponential rate, we see that P is correspondingly 
decaying at an exponential rate. Thus chaos (for Hamiltonian 
systems) necessarily implies the collapse of the phase space 
along certain directions. 



B. The Symplectic Eigenskeleton 



2) We define the matrices S and N as follows: 



In this section we expose a special basis associated with 
any hnear(ized) symplectomorphism that resists collapse. We 
shall refer to the symplectomorphism as <^t, keeping the 
dynamical setting (i.e. (I)t is the phase flow of a Hamiltonian 
system) in mind. 

Theorem 2 (The Symplectic Eigenskeleton): Consider a 
symplectomorphism (f>t : M ^ M that takes the initial 
point xo to (ptixo) = X. Let $ : T^gM T^M be the State 
Transition Matrix (STM) of the mapping. Let ^ = ■ ^ 
and let {^i, 771, ... , ^jv, iin} be the orthonormal eigenbasis 
of ^. Then the following are true: 

1) There is an interdependency amongst the vectors of 
The eigenvectors occur in pairs, where the can 
be taken to be 

Vi = J-^i 

where the associated eigenvalue of rji is A^^ if A, is 
the eigenvalue associated with ^j. 

2) The linear transformation T that takes the standard 
basis to the eigenbasis of 

T: {pi,qi,...,PN,qN} ^ {^i,Vi, ■ ■ ■ ,^n,Vn}, 

is symplectic. Moreover, the couples make 
symplectic pairs. 

3) The vectors {$ • ^1, $ • ryi, . . . , $ • ^n, $ • Vn} are 
orthogonal. Moreover, 



and ||*-?7,|| = V^Ap^ 



4) If a 2fc-dimensional symplectic subvolume T is ini- 
tially parallel to k of the eigenskelton planes, then the 
linearized transformation <I> preserves its volume, i.e 

Vol2fe(T) = Vol2fe($(T)). 

We call the symplectic eigenbasis of the matrix ^I* — $^ • $ 
the symplectic eigenskeleton of the transformation (j)t, as 
it is a property structure of the transformation which resists 
collapse over a discrete time t. 
Proof: 

1) Consider the i-th eigenvector of ^ with eigenvalue 
Xf. 

* • = A«e« 

Taking the transpose of this equation, right-multiplying 
by J - and then recognizing the identity ^'^ • J - = 
J, we see that 

Taking the transpose once more and multiplying by —1 
(whilst noting = ^ and = — J) we have 

Hence, the vectors r]i = J ■ £,i are also eigenvectors of 
with eigenvalues A^^. 



and N : 



where N = J ■ 'E from Part 1. The transformation 
matrix T can be represented as: 



T = 



We have temporarily reordered our representation of 
the basis, so that 



J 



Noting again that N — J ■ one easily sees 



Due to the orthonormality of the eigenvectors (i.e. • 
N = ■E = On sndE'^ ■E = ■ N = In), this 
expression reduces to: 

■J-T=J 

and hence the matrix T is symplectic. 
3) Renaming the eigenvectors of as {fijf^i ~ 
{^ijVi}iLi, we have from the orthonormality of the 
eigenbasis of 



Sin 



But 



($ • Vi) ■ • Vj) =vf Vj = XjSij 

Hence the vectors {$ • Vi}f^-^ are also orthogonal. 
Moreover: 



1$ 



Similarly 



1 



4) If T is a 2-dimensional area element spanned by ^1 
and rji, then the area expansion is 



Vol2($(T)) 
Vol2(T) 



1 



But since this is true of any area element initially paral- 
lel to one of the symplectic eigenskeleton planes, and 
the symplectic eigenskeleton 2k volumes are simply 
direct sums of these area elements, the result follows. 



VI. Applications to Control 

We view this paper as a theoretical paper which studies 
some of the fundamental constraints in the propagation of 
volumes rather than trajectories in dynamical and control 
systems. This idea has already been advocated in viability 
theory and in some robust control design papers, see e.g. 
Mayne [6]. In the final subsection below we present some 
future directions for using the theory presented here in the 
context of control. The examples below, though they do 
not utilize the full breadth of the theoretical developments 
presented in the paper, were chosen to illustrate some key 
ideas regarding propagation of surfaces and uncertainties in 
the control theory setting. 

A. The Kinematic Heisenberg System 

The Heisenberg System is a classical underactuated kine- 
matic control problem with nonholonomic constraints, see 
Bloch [3], Brockett [4]. The configuration manifold for the 
system is given by Q = R'^, with coordinates q = {x, y, z). 
Motion is constrained by the relation i = yx~xy. Supposing 
we have controls u, v over the x and y velocities, the 
kinematic control system can be written: 



X = u 

y ^ V 

z — yu — XV 



(11) 



Suppose the initial state of the system is given to be 
within the two-dimensional uncertainty distribution (t(0) — 
{{x,y,z) £9? : X e [-1,1], y G [-1,1], z = 0}, and we 
wish to determine an open loop control law that maneuvers 
the particle to the point (0, 0, 1) during the time interval 
t e [0, 1] in some optimal sense. 

Let a{t) be the time evolution of (t(0) to time 
t. We can parameterize the surface a{t) by the 
initial data {X,Y) e [-1,1] x [-1,1], so that, 
at time t, the surface is given parametrically by 
{x{X,Y;t),y{X,Y;t),z{X,Y;t)). The distance from 
an arbitrary point on the final surface cr(l) to the target 
point is ^x{X, F; 1)2 + y{X, Y; 1)2 + (1 - z{X, Y; 1))2. 
The dynamics ( fTTT i depend upon the choice of control 
(u(t), v{t)). We thus pose the following control problem: 

Problem: Choose {u{t),v(t)) so as to minimize: 



^(1) 



y 



(1 



(12) 



Solution: First we need to compute the determinant of the 
surface metric. The State Transition Matrix (STM) dynamics 
are given by: 






u 



■ $, $(0) = 







which can immediately be integrated to yield: 

m = 



-(2/(0- y(0)) (x(t)-x(O)) 



The metric determinant of the surface at time t is given by 
the Gram determinant of the first two columns of ^{t): 



g{X, Y;t) = l + ix{t) - Xf + {y{t) - Yf 
Thus (fT2l l is equivalent to: 



(13) 



^/x^ + y2 + (1 _ z)^y^g{X, Y; 1) dX dY 



1-1 J -1 
where {x,y,z) 



: {x{X,Y-l),y{X,Y-l),z{X,Y-l)). 
Without loss of generality, let us instead minimize: 

/ = j\ f ^ Hi? + y{l? + (1 - ^(1))'] .9(1) dXdY 

(14) 

Define now: 

— u{t) dr, 1^ — !'(''") dr, ct — [vu — fiv) dr 
Jo Jo Jo 

so that the solution to (fTTI) can be expressed as: 

x{t) = X + ii{t) 
y{t) = Y + y{t) 
z{t) = Y ^l{t) - Xv{t) + a{t) 

This exposes the dependence of x, y, and z on the initial 
conditions X and Y. Inserting into the surface metric ( fTSl ). 
one can explicitly integrate ( fT4b to find: 

f = ^{l + H^ + u^) (4^2 ^ 4j,2 _^ „ 6a + 5) 

This function has a global minimum at/i = 0, u — 0, a — I. 
Any control law that satisfies: 



u{t)dt = 0, / v{t)dt = 0, / (v{t)u(t)-^i{t)v{t))dt = 1 



will leave the final uncertainty distribution as close to the tar- 
get point as possible, in the above sense. Notice that a phys- 
ical interpretation of the vector quantity (/i(i), v{t), a{t)) is 
that it is the position vector of the point on the surface that 
was initially at (0, 0, 0). Thus, any control law that leaves the 
centroid of the surface at the target point will automatically 
minimize (fT2] l. One such trajectory, given in Bloch [3] using 
2/(0) = 0, is: 



Ki) = 
ait) = 



1 



/27r 
1 



sin(27rt) 
(1 - cos(27ri)) 



V2n 
t{l - sin(27rt)) 



The uncertainty surface at various time snapshots for the 
control law u{t) — fJ,{t), v{t) — v{t) is given in Fig. |2] 

B. The Falling, Rolling Disc 

Consider the falling rolling disc of radius r = 1, Fig. 
[51 whose configuration is described by the contact point 
{x^y) and the Classical Euler angles {(j),9,ip). Suppose we 
have direct control over the body-axis angular velocities 
u = (psinO, V = 0, w = (j)C0s6+ip, and suppose the system 



Since the STM ^{t) is initially the identity, we find: 



t=0.75 




2 -2 



Fig. 2. Optimal Uncertainty Maneuver 




Fig. 3. Euler Angles of the Falling Rolling Disc 



is subject to nonholonomic constraints x + ip cos <p ~ and 
y + ip sin (f) = 0. The dynamics is given by the system: 



X 

y 



u cot 9 cos 4> — w cos ( 
u cot 9 sin (f) — w sin 4 
ucsc9 

V 

—u cot + w 



Using the notation q — f{q,u), the coefficient matrix in the 
STM dynamics equation is: 



Of 



— (m cot 9 ~ w) sin </> 

(w cot 9 — w) cos (/) 









u csc^ 61 cos (f) 

M csc^ 9 sin (/) 

— ucot0csc6' 



—ucsc^9 



I A C 

1 B D 

1 £; 

1 

1 



where 

A = 

B = 

C = 

D = 

E = 

F = 



- {u cot 9 — w) sin dt 



{u cot 9 ~ w) cos (j) dt 



(^—{ucot9 — w)sm(f)E + ucsc'^6cos4)) dt 
(^{u cot 9 — w) cos (fiE + u csc^ 9 sin (j)^ dt 



-u cot 9 CSC 9 dt 



csc^ 9 dt 



Supposing there is initial uncertainty in the 9 and (j) compo- 
nents, it is the third and fourth columns of the STM that will 
be crucial in determining the uncertainty evolution. Suppose 
further our desire is that the projection of the final uncertainty 
onto the x, y planes has zero area; i.e. at worst there is a one- 
dimensional uncertainty in the contact point position. Then 
we wish at time t = 1, that A(l)i:>(l) - S(1)C(1) = 0. To 
achieve this, one may use any control law with that satisfies 
the relation 

u cot 9 — w = 

Such a control law will leave A{t) = 0, B{t) = 0, so that 
the uncertainty projection onto the x — y plane has zero area 
for all time. 

C. Future Directions in Control 

As the same symplectic constraints apply to the evolution 
of an optimal control system's states and co-states, these 
results also have an implication for the stability and robust- 
ness of an optimal feedback control law. This aspect of the 
study can be reduced to two fundamental approaches, the 
implication of initial value distributions on the subsequent 
evolution of a trajectory in the neighborhood of the true 
optimal trajectory, and how the symplectic invariants mani- 
fest themselves in the solutions of two-point boundary value 
problems. 

First, how do uncertainties in the initial state or in the 
initial application of the control map to the target conditions? 
As is well known, by definition an explicit optimal feedback 
control law is asymptotically stable when restricted to the 
state variables. However, as the necessary conditions from 
which the feedback control law can arise form a Hamiltonian 
system, this implies that the co-states are unstable and should 
diverge. This becomes an issue if the state is not perfectly 
determined or if the control function is not exact but only 



lies in a neighborhood of the true optimal control, and 
should lead to instabilities arising in the state variables of 
the system. These relationships can be studied using integral 
invariants and symplectic capacities to determine the robust- 
ness of the specific optimal control laws by studying how the 
phase volume surrounding them maps under the necessary 
conditions. Of special interest will be the identification of 
the maximum and minimum uncertainty growth directions. 

Second, given an optimal control feedback law (i.e., given 
the solution to the Hamilton-Jacobi-Bellman equation), how 
do simultaneous uncertainties in both the initial state and 
target state affect the distribution of the adjoints, and what 
structure may lie within these distributions that arise from 
the Hamiltonian formulation of the necessary conditions? 
Applying the Hamilton Principle Function approach, which 
provides an explicit solution to the two-point boundary 
value problem and which is directly, analytically related to 
the optimal control, we implicitly define an initial set of 
optimal controls that will lead to a proscribed final region in 
the neighborhood of the nominal target state. This defines 
for us an open set of controls, within which lie optimal 
trajectories that all achieve the final state to within some 
desired, and proscribed, accuracy. Depending on the size and 
the distribution of the initial uncertainties and the tolerable 
final uncertainties, we can identify a symplectic width which 
should provide expUcit ranges in the set of initial controls 
that will lead to a guaranteed, optimal arrival in the vicinity 
of the final state. Such a development can provide additional 
insight into the robustness of optimal controls and how 
gracefully they will degrade when we allow for finite miss 
distances for the target state. 

There is also a clear identification between distributions 
in phase space and probabilistic interpretations of the state 
of a system. Thus, our research also has a direct bearing 
on predicted uncertainties in a dynamical system after being 
mapped in time, and will define for us an absolute minimum 
region within which the uncertainty of the system can be 
isolated. 

VII. Conclusion 

We showed how the expansion of subvolumes in the local 
neighborhood of a nominal trajectory leads to the local col- 
lapse of the supporting phase space. Moreover, we produced 
a preferred basis, the symplectic eigenskeleton, which resists 
collapse and returns uncertainty distributions that are initially 
parallel to the basis to their minimal uncertainty state at a 
fixed final time. 
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